Conformal geometry of the (discrete) 
Schwarzian Davey-Stewartson II hierarchy 



B.G. KONOPELCHENKO 

Dipartimento di Fisica, Universita di Lecce and Sezione INFN, 73100 Lecce, Italy 

W.K. SCHIEF 

School of Mathematics, The University of New South Wales, 
Sydney, NSW 2052, Australia 

Abstract 

The conformal geometry of the Schwarzian Davey-Stewartson II hier- 
archy and its discrete analogue is investigated. Connections with discrete 
and continuous isothermic surfaces and generalised Clifford configurations 
are recorded. An interpretation of the Schwarzian Davey-Stewartson II 
flows as integrable deformations of conformally immersed surfaces is given. 

1 Introduction 

Due to the (re-) discovery of a variety of important connections between the 
differential geometry of surfaces and integrable systems, (classical and modern) 
differential geometry has been widely recognised as an integral part of soliton 
theory (see, e.g., [l]-[3]). However, the fundamental nature of geometry in 
the context of integrable systems is a subject of ongoing research and recent 
investigations have uncovered unexpected geometric links. For instance, it has 
been established that Hirota's master equation [4] in its Schwarzian form and 
the associated scalar Schwarzian Kadomtsev-Petviashvili (SKP) hierarchy are 
encapsulated in Menelaus' classical theorem of plane geometry [5] -[7]. 

In the present paper, we embark on a study of the geometry of the Schwarzian 
Davey-Stewartson II hierarchy and its discrete analogue, the quaternionic dis- 
crete SKP (qdSKP) equation. We demonstrate that the qdSKP equation and 
various associated continuum limits are canonical objects of conformal (Mobius) 
geometry in 1R 4 . In particular, we establish important connections with both 
discrete and continuous isothermic surfaces and generalised Clifford point-circle 
configurations. We also show that the Schwarzian Davey-Stewartson II hierar- 
chy explicitly defines integrable deformations of conformal immersions in R 4 . 

2 The multicomponent discrete Schwarzian KP 
equation 

The multicomponent KP hierarchy houses a variety of important soliton equa- 
tions such as the Davey-Stewartson and N-w&ve equations and their associated 



hierarchies [8]. The Schwarzian KP (SKP) hierarchy consisting of the singular- 
ity manifold equations for the multicomponent KP hierarchy has been shown to 
admit an elegant compact formulation [9]. Indeed, it has been established that 
if an N x N matrix 

$(t), t=(t 1 ,t 2 ,t 3 ,...) (1) 

depending on an infinite number of 'times' t n constitutes a solution of the SKP 
hierarchy then the six solutions 

$i=Ti$, $ jfc =T,T fe $, i,k= 1,2,3, i ? k, (2) 

where the 'shift' operators Tj are defined by 

m(t) = S(t + [Oj]) = + Oi, t 2 + y . *3 + • ■ •) (3) 

and aj = const, obey the algebraic 6-point relation 

M($i,$i2,$2,$23,*3,*3l) = -1- (4) 

Here, the multi-ratio M of six matrices P 1 , . . . , P e is defined by 
M(P\P 2 ,P 3 ,P 4 ,P 5 ,P 6 ) 

= (P 1 - P 2 )(P 2 - P 3 )" 1 ^ 3 - P 4 )(P 4 - P 5 )- 1 ^ 5 - P 6 )(P 6 - P 1 )- 1 . 

(5) 

The entire SKP hierarchy may then be retrieved from (4) by considering a 
canonical limit in which a, — > [9]- [11]. 

As indicated above, the multi-ratio relation (4) represents an algebraic su- 
perposition formula for six solutions of the multicomponent SKP hierarchy. As 
in the scalar case [5, 10], it also constitutes an algebraic relation ('permutabil- 
ity theorem') for six solutions generated by a variant of the classical Darboux 
transformation [12, 13]. Iterative application of the Darboux transformation 
then produces lattices of solutions of the multicomponent SKP hierarchy. In 
particular, the multi-ratio relation (4) may be interpreted as an equation de- 
fined on a H? lattice. Thus, in the following, we regard (4) as a discrete equation 
for a matrix-valued function 

$ : 1? -> C^, (m, n 2 , m) ^ n 2 , n 3 ), (6) 

where the indices on $ denote translations on the lattice, that is, for instance, 

$ = $(ni,n 2 ,n 3 ), $i = $(ni + 1, n 2 , n 3 ), $ 23 = $(m, n 2 + l,n 3 + 1). (7) 

The discrete equation (4) has come to be known as the multicomponent 
discrete SKP (dSKP) equation since it encodes the complete multicomponent 
SKP hierarchy. Indeed, the genesis of the multicomponent dSKP equation im- 
plies that any member of the multicomponent SKP hierarchy may be obtained 
from (4) by applying an appropriate continuum limit. The integrable nature of 
the multicomponent dSKP equation is inherited from both the original deriva- 
tion (l)-(4) and the construction via Darboux transformations. The following 
analysis is concerned with the conformal geometry of the quaternionic dSKP 
equation and its various continuum limits. 
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3 The quaternionic dSKP equation 



An important property of the multicomponent dSKP equation is its invariance 
under fractional linear transformations of the form 

$ -y = (A® + B){C<f> + D)- 1 , (8) 

where A, B, C and D are arbitrary constant matrices. It is therefore natural to 
investigate whether contact may be made with conformal (differential) geometry. 
Thus, in the remainder of the paper, we assume that the matrix $ takes values 
in the space of quaternions H. The latter is identified with a four-dimensional 
Euclidean space R 4 via 

R 4 3 (a, b, c, d) <-» {at + bi + cj + dk) E JH, (9) 

where the matrices 1, i, j, k are defined by 

'-(;:)• >=(->»)' j-tf-o 1 )- -(o 1 ?)- <«•) 

In particular, this isomorphism gives rise to the identities 

X 2 =dctX, XX* = (detX)l, X Y = ^tv{XY^) (11) 

for any quaternions X, Y E H and their vectorial analogues X,Y E R 4 . 
Accordingly, in the case $ E H, we may refer to (4) as the quaternionic 
dSKP (qdSKP) equation. The invariance (8) encodes the group of orientation- 
preserving conformal transformations in R 4 provided that A, B, C and D are 
quaternions. 

3.1 The 'geometric' continuum limit 

As alluded to in the preceding section, the complete quaternionic SKP hierarchy 
may be retrieved from the qdSKP equation via appropriate sophisticated limits. 
We embark on a study of the geometric implications of this fact in the next 
section. Here, by contrast, we focus on the natural 'geometric' limit in which 
the differences Ai<I> = $^ — $ are regarded as approximations of derivatives, 
that is 

$ l = $ + e$^+0(e 2 ), i= 1,2,3, (12) 

where e is a lattice parameter and $ Xi = d^/dxi. In the limit e — > and 
{xi, X2, £3) = {x,y,z), the qdSKP equation reduces to 

= 1. (13) 

The latter becomes an identity if the symmetry 

% + = (14) 
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is imposed. The scalar analogue of this constraint has been shown to lead to 
conformal maps or (anti-)analytic functions on the complex plane [5]. As in the 
scalar case, in order to obtain a non-trivial continuum limit from the qdSKP 
equation subject to the constraint (14), it is required to take into account second- 
order terms in the Taylor expansion (12) and consider the terms in the expansion 
of the qdSKP equation which are linear in e, that is 



* x z y x z 



*yz 



y x z 



(15) 



1 



<&xz + ^xx - 1 \$xz + 77$ 



sr 1 =0. 



Simplification by means of (14) then yields 



(16) 



By virtue of the isomorphism R 4 = H, any solution of the quaternionic equa- 
tion (16) gives rise to a surface immersed in 1R 4 and parametrised by the co- 
ordinates x and y. Its position vector r = (r, r) G R 4 is obtained from the 
decomposition 

$ = rl + r-e, (17) 

where the 'vector' e is defined by e = (i,j,k). In terms of r, equation (16) is 
readily seen to translate into 



_ • y • yy 



Vx ' r yy r i Vx ' r y r 



yy 



(18) 



Thus, in the natural geometric continuum limit, the qdSKP equation subject 
to the symmetry constraint (14) governs surfaces the position vector of which 
obeys the second-order equation (18). 

By construction, the surfaces defined above belong to conformal differential 
geometry (see, e.g., [14]) due to the invariance of (18) under the group of con- 
formal (Mobius) transformations as induced by (8). These are integrable in the 
sense that they are given in terms of solutions of the soliton equation (16). In 
order to proceed, it is now observed that (16) implies that 



Qy = \[^yy^y\Qi 



where the quaternion Q is defined by 

It may therefore be admissible to impose the constraint 



(19) 

(20) 
(21) 



4 



since, in this case, relation (19) is identically satisfied. In the generic case, that 
is 76 <f> y , the constraint (21) is equivalent to the pair 

tr($ K $- 1 )=0, det^*" 1 ) = 1 (22) 

which, in terms of the position vector r, becomes 

r 2 x =r 2 y , r x -r y = 0. (23) 

The latter conditions may be used to cast (18) into the form 

r x ■ r xy r y ■ r xy 

I xy - 2 V x -f T y [Z'i) 

x y 

or, equivalcntly, 

r xy = ar x + br y , (25) 

where the real coefficients a and b arc determined by the constraints (23). Thus, 
it has been established that the surfaces defined by the system (23), (25) con- 
stitute a subclass of the surfaces associated with the continuum limit of the 
qdSKP equation considered here. 

The geometry of the reduction (23), (25) is now readily revealed. Firstly, the 
constraints (23) are equivalent to demanding that the first fundamental form of 
a surface be conformally flat, that is 

dr 2 = n(dx 2 + dy 2 ), n = r 2 x = r 2 y . (26) 

The coordinates x and y are therefore conformal coordinates [15]. Secondly, 
the hyperbolic equation (25) expresses the fact that the vector r xy is tangential 
to the surface so that the coordinates x and y are conjugate [15]. A surface 
which may be parametrised simultaneously in terms of conformal and conjugate 
coordinates is termed an isothermic surface.' 1 ' Hence, we conclude that the 
continuum limit (16) of the qdSKP equation subject to the constraint (21) is 
associated with isothermic surfaces. The latter are classical and have been 
investigated extensively with respect to both geometry and integrability (see, 
e.g., [2, 3] and references therein). The connection with isothermic surfaces 
therefore provides a first indication of the fundamental nature of the qdSKP 
equation. 

3.2 The conformal geometry of the qdSKP equation. Dis- 
crete isothermic surfaces 

The qdSKP equation has recently been given a geometric interpretation [16] in 
terms of a novel generalisation of Clifford's classical C4 point-circle configuration. 
A C4 Clifford configuration is constructed in the following manner [17, 18]: 
Consider a point P° on a plane and four generic coplanar circles S 1 , S 2 , S 3 , S* 4 
passing through P°. The additional six points of intersection are labelled by 
P 12 , P 13 , P 14 , P 23 , P 24 , P 34 , where the indices on P lk correspond to those of 
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Figure 1: A C4 Clifford configuration 



the intersecting circles S z and S (cf. Figure 1). Any three circles S z ,S k ,S l 
intersect in three points and therefore define a circle S lkl passing through the 
points of intersection p* k p tl ^p kl Clifford's circle theorem then states that, 
remarkably, the four circles S 123 , S 124 , S 134 , S 234 meet at a point p 1234 . It is 
noted in passing that Clifford configurations (C n ) exist for any number of initial 
circles S , . . . , S n passing through a point P . 

In [5] , it has been shown that the six points P %k of a C4 Clifford configuration 
are algebraically related by the multi-ratio relation 

M(P 14 ,P 12 ,P 24 ,P 23 ,P 34 ,P 13 ) = -1 (27) 

if the plane is identified with the complex plane so that the points P tk are 
regarded as complex numbers. The above multi-ratio relation encodes nothing 
but Mcnelaus' classical theorem of plane geometry if the point P° is mapped 
to infinity by a conformal transformation [19]. Here, it is important to note 
that if we set aside the points P° and p 1234 then a C4 Clifford configuration 
exhibits the combinatorics of an octahedron if the six points P lk and eight circles 
S l , S lkl are identified with the vertices and faces of an octahedron respectively. 
In particular, the multi-ratio relation (27) admits the full symmetry group of 
an octahedron acting on the entries P tk . Moreover, Zicgenbein [18] has proven 
that all circles and points of a C4 Clifford configuration appear on equal footing 
in the sense that the angles made by four oriented circles passing through a 

t Accordingly, a surface is isothermic if and only if its curvature coordinates arc conformal 
modulo a suitable rcparamctrisation x — > f(x), y — > g(y) [15]. 
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point are the same for all eight points. The converse of Ziegenbein's theorem is 
also valid. 

It has been established in [5] that any 'generic' six points P lk on the complex 
plane belong to a C4 Clifford configuration if and only if the multi-ratio relation 
(27) is satisfied. It is therefore natural to inquire as to the geometric significance 
of the quaternionic multi-ratio condition 

M(P 14 , P 12 , P 24 , P 23 , P 34 , P 13 ) = -1, (28) 

where the quaternions P %k are regarded as points in a four-dimensional Eu- 
clidean space. It is evident that for any 'generic' six points P %k G K 4 , the 
circles S l and S tkl may still be constructed but, in general, these do not in- 
tersect in any points P° and p 1234 . However, if one imposes the Ziegenbein 
property then the six points P lk can no longer be arbitrary. In fact, it turns 
out [16] that the Ziegenbein property is equivalent to the quaternionic multi- 
ratio condition. t Thus, the latter condition gives rise to generalised C4 Clifford 
configurations if their definition is based on the Ziegenbein property. 

The preceding discussion of Clifford configurations implies that the qdSKP 
equation 

M($i, $12, $ 2 , $ 23 , $3, $31) = -1. (29) 

enshrines 'collections' of generalised C4 Clifford configurations. This has been 
made precise in [16] . On the other hand, the connection with classical isothermic 
surfaces as established in Section 3.1 raises the question as to whether the stan- 
dard integrable discretisation of isothermic surfaces [2] is related to generalised 
Clifford configurations. A quadrilateral lattice (discrete surface) 

$ : 1? JR 4 ^ JH (30) 

is termed isothermic if the quaternionic cross-ratio 

Q($, $1, $12, $2) = (* - $i)(*i - *i2) _1 (*i2 - $ 2 )($ 2 - (31) 

associated with any quadrilateral obeys 

Q($,$i,$i2,$ 2 ) = -1. (32) 

Since the above quaternionic cross-ratio condition is invariant under conformal 
transformations and any four points in R 4 may be mapped to the plane by 
means of an appropriate conformal transformation, the quadrilaterals are in- 
scribed in circles and their classical (scalar) cross-ratio is —1 [2]. This is the 
geometric content of the quaternionic cross-ratio condition (32). It is observed 
that the quaternionic cross-ratio condition constitutes a natural discretisation 
of the constraint (21). However, since the quadrilaterals are planar, it also 
discretises the hyperbolic equation (25) so that (32) may indeed be regarded 
as a discrete version of the conditions (23), (25) defining classical isothermic 
surfaces [2]. 

t Modulo an inversion with respect to a hypersphere. 
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The canonical discrete analogue of the constraint (14) is given by the trans- 
lational symmetry 

$ 23 = $. (33) 

The latter may be used to eliminate quantities which carry an index 3 from the 
qdSKP equation (29). On rearranging terms, one obtains 

Q(*i, *12, *2, *) = T- 2 Q(® 12 , $1, *, $ 2 ), (34) 

where T 2 f{n-i) = f(n 2 — 1). If the quaternionic cross-ratio condition (32) holds 
then both cross-ratios in (34) are —1 due to the symmetry group of the cross- 
ratio condition. Accordingly, the pair (32), (33) constitutes a reduction of the 
qdSKP equation. In geometric terms, this implies that any discrete isothermic 
surface extended to a three-dimensional lattice via the translational symmetry 
(33) represents a collection of (degenerate) generalised C4 Clifford configurations 
as introduced in [16] and alluded to in the preceding. 

4 Deformation of conformal immersions induced 
by the Schwarzian Davey-Stewartson II hier- 
archy 

In order to investigate the conformal differential geometry of the continuous 
qSKP hierarchy, we now recall the connection between the qdSKP equation 
and the (adjoint) eigenfunctions associated with the quaternionic discrete KP 
equation, that is the discrete Davey-Stewartson II equation [9]-[ll]. Thus, it is 
readily verified that the linear system 

02 = ^i(A 1 $)- 1 A 2 $, <f> 3 = <^ 1 (A 1 $)- 1 A 3 $ (35) 

for a quaternionic function <f> is compatible modulo the qdSKP equation (29). 
Hence, one may introduce a quaternionic function of 'separation' ip according to 

A,$ = (36) 

The compatibility conditions [A i; Aj,]$ = then yield 

i4> i -<t>k)4>* 1 = ii>- 1 y>k-ii>i). (37) 

On the one hand, addition of the three relations (37) produces 

(</>! - ^fe 1 + (02 ~ 03)023 + (03 ~ 0l)031 1 = ( 38 ) 

which constitutes the eigenfunction equation for the quaternionic discrete KP 
equation. The latter represents nothing but a discrete Davey-Stewartson II 
equation so that (38) may be regarded as a discrete version of the modified 
Davey-Stewartson II equation which is known as the Ishimori equation [20] . On 



8 



the other hand, elimination of the cigcnfunction <j> from (37) in a similar manner 
gives rise to the 'adjoint' eigenfunction equation 

V'l" 1 ^! - Vl2) + ^ 2 " 1 (V'12 - V>2 3 ) + ^(fe - ^3l) = (39) 

with i/> being an adjoint eigenfunction of the discrete Davey-Stewartson II equa- 
tion. It is emphasised that the qdSKP and (adjoint) quaternionic discrete KP 
eigenfunction equations are equivalent. For instance, if (f> is an eigenfunction 
obeying (38) then the linear system (37) for ip is compatible and ip constitutes a 
solution of the adjoint eigenfunction equation (39). Moreover, by construction, 
the existence of a function $ satisfying the defining relations (36) is guaranteed 
and $ is indeed a solution of the qdSKP equation (29). 

The continuum limit to the qSKP hierarchy requires the introduction of the 
gauge transformations 

4>^A^A^A^^ ^^^V 1 , (40) 

where the constant matrices Ai constitute non-degenerate diagonal quaternions. 
In the limit in which the qdSKP equation reduces to the nth-order qSKP equa- 
tion, the system (36) becomes [9] -[11] 

* x = ^Mi0, § y = ipA 2 4>, $ t =B„ftM), (41) 

where B n is bilinear in ip, <f> and their derivatives. If we make the choice A± = 1 
and A 2 = k and eliminate the eigenfunction <p from (41) then we obtain the 
pair 

$ y = N$ x , $ t = £,($), (42) 

where 

N = i/'kV'" 1 (43) 

and £„($) = B n (ip, ip~ 1 $ x ). Accordingly, C n constitutes a linear operator 
acting on $ with coefficients depending on ip and its derivatives. It turns out 
that these coefficients may be written as differential expressions in N and a 
scalar auxiliary function ip [9]- [11]. The compatibility condition for the linear 
system (42) then gives rise to the nth-order Ishimori equation. For instance, in 
the case n = 2, the Lax pair (42) assumes the form 

$ y = N1> x , * t = N(& xv -<py$ x -<p x $ y ) (44) 

and its compatibility condition yields 

N t = NN X y + NyN X + lf X N X - IfiyNy - (ifi ^ + Ifiyy) N . (45) 

Now, the definition (43) of TV implies that 

triV = 0, N 2 = -l, N = S-e, S 2 = 1 (46) 
and hence decomposition of the matrix equation (45) produces 

S X S X y -\- (p x S X (fiySy 

(47) 

<Pxx + fyy + (S x X S y ) • S = 0. 
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The latter represents the Ishimori equation which was first set down in [20]. It 
is interesting to note that the 'topological charge' 

Q= i^/ (SxXSy)-Sdxdy (48) 

is preserved by the Ishimori flow (47)i. In fact, the topological charge may be 
shown to be invariant under all higher-order Ishimori flows. It is also remarked 
that the Ishimori hierarchy is amenable to the inverse spectral transform (1ST) 
method [21, 22]. 

The geometry of the Schwarzian Davey-Stewartson (SDS) II hierarchy en- 
coded in the Lax pair (42) for n = 2, 3, . . . is unveiled by focussing on the 
'scattering problem' (42) i (cf. [23]). Thus, if, for any fixed t, we identify the 
quaternionic function $ with the position vector r of a surface in R 4 then the 
properties (46)i,2 show that the metric of the surface is once again given by (26) 
with x and y being conformal coordinates. Any SDSII flow (42)2 therefore de- 
fines an intcgrable deformation of conformal immersions of surfaces in 1R 4 , where 
the independent variable t is regarded as the deformation parameter. The asso- 
ciated Ishimori flows possess an infinite set of invariants, the simplest of which 
is given by the topological charge (48). Moreover, the invariance of the qdSKP 
equation under the Mdbius transformation (8) guarantees that the SDSII de- 
formations are covariant under conformal transformations in R 4 . For instance, 
the Lax pair (44) is form-invariant under 

jV->$ _1 iV$, ^^^-lndet^ (49) 

corresponding to the composition of an inversion and a reflection in R 4 given 

by 

An important feature of the immersion (42) i and the deformations (42) 2 is 
that these are formulated explicitly in terms of geometric quantities, namely the 
position vector r (or <j>) and the quantity N . The latter is known as the 'left 
normal' of the surface [23]. If tr$ = 0, that is if the surface is embedded in 
R 3 = ImH, then the immersion formula (42) i decomposes into 

x v = S x r x , S • r x = (51) 

so that S constitutes the unit normal 

S = ^ (52) 

to the surface. However, in general, N is not normal to the surface and the 
constraint tr $ = is not preserved by the SDSII flows (42) 2 . In this connection, 
it is interesting to investigate the stationary points of the SDSII deformations. 
For instance, if $ t = then the flow (44)2 reduces to 

$ xy = <p y $ x + Lp x $ y (53) 
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which is nothing but the hyperbolic equation (24) with 

p=^lnfi. (54) 

Thus, the class of surfaces which is preserved by the deformation associated 
with the Schwarzian Davey-Stewartson II equation coincides with the class of 
isothermic surfaces in R 4 . The metric of any such surface is given by 

dr 2 = e 2ip (dx 2 + dy 2 ) (55) 

which affords an immediate geometric interpretation of the auxiliary function ip. 
Moreover, the auxiliary equation (47)2 constitutes the Gaufi equation [15] 

fxx + fyy + ICe 2v = 0, (56) 

where the Gaufiian curvature /C of the surface is given by 

IC = (S x xS y )-Se- 2v . (57) 

The topological charge Q defined by (48) is therefore proportional to the total 
Gaufiian curvature. Finally, it is noted that the constraint tr <f> = is compati- 
ble with (53) and gives rise to classical isothermic surfaces in R 3 . 
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